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We study the effect of disordered ripples on the conductivity of monolayer graphene flakes. We 
calculate the relaxation times and the Boltzmann conductivities associated with two mechanisms. 
First, we study the conductivity correction due to an external in-plane magnetic field B\\ . Due to 
the irregular local curvature found at graphene sheets deposited over a substrate, B\\ can be mapped 
into an effective random magnetic field perpendicular to the graphene surface. Second, we study the 
electron momentum relaxation due to intrinsic pseudo magnetic fields originated from deformations 
and strain. We find that the competition between these mechanisms gives rise to a strong anisotropy 
in the conductivity tensor. This result provides a new strategy to quantitatively infer the strength 
of pseudo-magnetic fields in rippled graphene flakes. 

PACS numbers: 72.80.Vp,73.23.-b,72.10.-d 


I. INTRODUCTION 

The electronic transport properties of bulk graphene 
are quite remarkable 1-4 . At room temperature graphene 
has mobilities 5 as high as fi ^ 15,000 cm 2 V -1 s -1 , that 
are significantly larger than any other semiconductor. At 
low temperature, however, the typical mobilities increase 
only up to /i ^ 200,000 cm 2 V -1 s -1 , which is a disap¬ 
pointing figure as compared with high quality GaAs het- 
erostructures. These observations triggered an intense 
theoretical and experimental activity to understand the 
disorder mechanisms that limit the mobility in graphene 
(see, for instance, Refs. 1-4 for a review). The theoretical 
understanding, acquired from the analysis of the Boltz¬ 
mann equation 1,3 , numerical simulations 2,6,7 , and field 
theoretical techniques 8 is that extrinsic disorder such as 
ad-atoms absorbed in the graphene surface, substrate 
charge inhomogeneities, and intrinsic disorder such as va¬ 
cancies play a key role. 

Motivated by recent experiments 9-12 , we study the ef¬ 
fect of extrinsic and intrinsic random gauge potential dis¬ 
order due to strain. We show that, although unlikely to 
be dominant in graphene deposited over standard sub¬ 
strates, these kinds of disorder give unique and sizable 
contributions to the conductivity. For high-quality sub¬ 
strates, there are even experimental evidences 12 that ran¬ 
dom strain can be a good candidate for the leading elec¬ 
tron relaxation mechanism in on-substrate graphene. 

Standard electronic transport experiments use samples 
where graphene flakes are deposited over an insulating 
substrate. In such setting, it has been experimentally es¬ 
tablished that the graphene surface is characterized by 
disordered static ripples 13-16 . For Si02 substrates, the 
latter have typical lengths of A ^ 5 — 30 nm and charac¬ 
teristic heights of h rrns ~ 0.2 — 0.5 nm. 

Lattice deformations due to ripples change the distance 
between the atoms in the graphene sheet. At the quan¬ 
tum level, site lattice displacements change the orbital 
bonding between the corresponding atoms, modifying the 
electronic structure of the material. In graphene, whose 


low-energy electronic properties are nicely described by 
a nearest-neighbor tight-binding model 1 , the occurrence 
of ripples change the tight-binding hopping terms 1,17,18 . 
For distortions with length scales much larger than the 
lattice parameter, characteristic to most samples 13-16 , 
it has been shown that the tight-binding model can 
be mapped into an effective Dirac Hamiltonian with a 
pseudo-magnetic vector potential 1,19 21 that depends on 
the lattice distortions. Hence, random ripples give rise 
to an intrinsic random gauge potential. The experimen¬ 
tal evidences of pseudo-magnetic fields are scarce and 
indirect, but quite remarkable. Strain fields have been 
invoked to associate the local density of states observed 
in graphene nanobubbles 22 to Landau levels with ener¬ 
gies corresponding to very high magnetic fields. To the 
best of our knowledge, no transport experiment has yet 
observed manifestations of this physical picture. 

Random magnetic fields can also be achieved by realiz¬ 
ing that disordered ripples in graphene and the roughness 
at semiconductor interface heterostructures share several 
common features. Starting at the late 80’s, a number of 
ingenious methods where used to characterize the inter¬ 
face roughness in the heterostructures 23-25 . One idea 
is particularly suited for graphene studies: By applying 
a strong magnetic field Hy, aligned to the plane of a 
heterostructure interface confining the two-dimensional 
electron gas (2DEG), the (smooth) interface roughness 
disorder gives origin to a local random magnetic field 
perpendicular to the 2DEG surface. The analysis of the 
electronic transport properties as a function of the ap¬ 
plied magnetic field gives quantitative information about 
the interface roughness. 

This setting was nicely explored by a recent 
experiment 9 , that combined information of the average 
conductivity and its weak localization correction 26,27 in 
graphene to extract the sample characteristic A and h rms . 
This procedure have been also used in the experimental 
study of other graphene systems 10,11 . Theory 28 shows 
that an applied B\\ on a rough surface gives rise to an 
effective dephasing and to the suppression of the weak 
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localization peak. In addition to this quantum correction, 
the random magnetic field due to B\\ also contributes to 
the electron momentum relaxation, which at high doping 
is accounted for by the Boltzmann theory 9 . This nice 
analysis does not consider the effect of intrinsic pseudo- 
magnetic fields due to strain, discussed above. 

Our focus is different. We study the combined effect 
of intrinsic and extrinsic random magnetic fields in the 
Drude conductivity. We revisit the analysis of the Boltz¬ 
mann equation in graphene 1 ,2,4,29 and calculate the con¬ 
tributions of random magnetic fields to the Drude con¬ 
ductivity. We show that the conductivity corrections due 
to an applied in-plane magnetic field B\\ on a rippled 
graphene flake depend on the direction of B\\ and are 
very anisotropic. We find that this result can be recon¬ 
ciled with the experimental findings of Ref. 9 by theo¬ 
retically treating the effect of strain and B\\ at the same 
footing. We also show that the combined effect of both 
sources of random magnetic fields provides a new exper¬ 
imental path to quantitatively probe the effects of strain 
fields in the low-energy electronic dynamics of graphene. 

The paper is structured as follows: In Section II we 
present the disorder model we employ to describe ripples 
and the corresponding effective Hamiltonian. We briefly 
explain the origin of an effective random field due to an 
applied in-plane magnetic field B\\, as well as from the 
intrinsic strain field due to ripples. In Sec. Ill we cal¬ 
culate the contributions to the Drude conductivity due 
to B\\ and strain using the Boltzmann transport equa¬ 
tion. We show how to account for the disorder potential 
anisotropy and discuss its consequences comparing with 
experiments. Finally, in Sec. IV we present our conclu¬ 
sions and an outlook. 


II. MODEL HAMILTONIAN 

In this section we present a model to study the effect 
of extrinsic and intrinsic sources of a random magnetic 
fields in the dynamics of electrons in corrugated graphene 
monolayer samples. 

Close to the charge neutrality point, the electronic dis¬ 
persion relation of pristine graphene monolayers is linear 
and has two degenerate components, with corresponding 
K and K' valley indices 1 . In the presence of a magnetic 
field, the effective electronic Hamiltonian for the A-valley 
reads 

H k = v F a • [p + eA(r)\ = H* + V(r), (1) 

where the vector potential A(r) has been included in 
H k by minimal coupling. Here v F « 10 6 m/s, a are the 
Pauli matrices acting on the sublattice space, and p is 
the electron momentum operator. The Hamiltonian for 
the K' valley has a similar structure 1 . 

In this description, a generic long-ranged disorder po¬ 
tential V(r) is represented in both K and K' valleys by 


where i = 0 stands for scalar disorder (with a o = I 2 ) 
and i = 1,2 for vector potential disorder, while i = 3 
represents a mass term. The focus of our study are (in¬ 
trinsic and extrinsic) disordered gauge fields, associated 
with V& and V< 2 ). In this paper we do not consider 
scalar disorder. 

Let us introduce 
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is an eigenstate of H/f, 0 = tan 1 (k y /k x ), s indicates 
particle (s = +1) or hole (s = — 1) doping, and 

V^ k , = jJdre^ k - k '> r V^(r) (5) 


is the momentum representation of Since we deal 

with elastic processes, we assume in the remaining of the 
paper that s = s '. 

At low temperatures, scalar disorder (short and long 
ranged) is the main source of momentum relaxation in 
graphene systems 2,4 . In this paper we use a phenomeno¬ 
logical transport time r s to account for effects of scalar 
disorder in the conductivity. We assume that r s is much 
shorter than the characteristic transport times due to 
random gauge fields. In Sec. IV, where we compare our 
results to experiments, we show that r s indeed dominates 
the conductivity in graphene, but some transport proper¬ 
ties are only explained by including effects due to random 
gauge fields. 

The ripple disorder model employed in this study is 
defined as follows: We describe the graphene sheet sur¬ 
face by 2 = h(r), where h is the surface displacement 
with respect to the reference plane 2 : = 0 at the posi¬ 
tion r = (x,y). The average of h is set to zero. In 
line with the experiments on graphene deposited over 
a substrate 9,13-16 , we further assume that the typical 
heights h rms are much smaller that the ripple lengths A. 

We model the ripple fluctuations in h(r) by the corre¬ 
lation function 


{h(r)h(r')) = h? rms F 


44 


(6) 


where (• • •) denotes an average over disorder. Al¬ 
though theory predicts a power law height-height correla¬ 
tion function for free-standing membranes 30 , experiments 
support single-parameter correlations for the (static) rip¬ 
ples of graphene deposited over a substrate. For latter 
convenience, let us define 


V(r) = 
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where A is the sample size. In reciprocal space 

(h(q)h{q')) = h 2 ms F(q)5 q _ q <, (8) 

where F(q) is the Fourier transform of the correlation 
function F(\r — r'\). 

We address two mechanisms that generate random 
magnetic fields. First, we study the case of an external 
strong magnetic field B\\ applied parallel to the graphene 
sheet. We show that, due to the ripples, B\\ gives rise to 
a random effective magnetic field F> ext (r) perpendicular 
to the graphene surface. Next, we discuss the intrinsic 
pseudo-magnetic field B- int originated by the strain field 
corresponding to the graphene sheet profile height h(r). 

A. Random magnetic field due to ripples and an 
in-plane external B-field 


and is expressed, in a convenient gauge for B\\ = By*, 
by the vector potential 

A x (r) = 0 and A y (r) =—B\\h{r). (13) 

Figure 2(a) illustrates a typical disorder realization 
of h(r) with fluctuations characterized by the Gaussian 
correlation function F(x) = exp (— x 2 /2X 2 ). The corre¬ 
sponding magnetic field B ext (r), normal to the graphene 
sheet, is shown in Fig. 2(b). While h(r) displays an 
isotropic disorder, R ext (r) is clearly anisotropic. The 
anisotropy direction of B ext (r) depends on the orienta¬ 
tion of B \\. 

The anisotropy is quantified by inspecting the autocor¬ 
relation function 

(14) 


Let us first consider the setup of a magnetic field ap¬ 
plied parallel to the sample z = 0, that has been exper¬ 
imentally investigated in a variety of systems 9,23-25 . For 
notational convenience, in what follows we fix the direc¬ 
tion of B\\ along the x-axis, namely, B\\ = B\\x. 


z 



FIG. 1. Sketch of h(r) along the x direction. The ripple 
amplitudes 5h are enhanced and made comparable with A to 
help the illustration. For convenience, we take B n = B\\x. 


As illustrated in Fig. 1, the parallel magnetic field B\\ 
has a component perpendicular to the surface 2 = h(r) 
that is given by 


B ext (r) =-B v n(r). (9) 

At the point Vq = ( xo,yo ), the surface z = h(r) has a 
unit normal vector 



We assume that the typical displacement magnitude is 
characterized by Sh. For A, we write 


71 ( 7 * 0 ) ~ (dh/dx,dh/dy,— 1) 


r=r 0 


(ii) 


Hence, the effective local perpendicular magnetic field 
reads 


Bextir) = —B\\ ■ Vh(r), (12) 


that can be expressed in terms of F by direct differentia¬ 
tion. Alternatively, going to reciprocal space, one writes 


dh{r ) dh(r') 
dx dx' 


= h 2 ims J2QlF(q)e- iq<r - r,) 


d 2 

— —h Ynis F(p), 


(15) 


with p = r — r '. 

Let us calculate (R ext (r)R ext (r / )) for the case where 
h(r) is characterized by a Gaussian correlation function, 
the same as in Fig. 2(a). The corresponding R ext (r) au¬ 
tocorrelation function reads 


(R ext (r)R ext (r')) 



A 2 



(16) 

where a is the angle between B\\ (or the x-axis) and p. 

Figure 3(a) shows the R ext (r) autocorrelation func¬ 
tion obtained by averaging over 10 5 ripple disorder real¬ 
izations of h(r), as defined by Eq. (6) with a Gaussian 
correlation function F (for more details see, for instance, 
Ref. 6). As expected, it coincides with Eq. (16) and ex¬ 
presses the anisotropy captured by a visual inspection of 
Fig. 2(b). 


B. Pseudomagnetic field due to strain 


The out of plane deformations of a rippled membrane 
described by h(r) can be associated with the strain tensor 
Uij(r ) given by 19,31,32 
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For simplicity we have neglected the effect of in-plane 
deformations. 
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FIG. 2. (color online) Typical disorder realization of (a) h(r), 
characterized by a Gaussian correlation function F, the cor¬ 
responding (b) B ex t(r) for an external B\\ applied along the 
x direction, defined in Eq. (12), and (c) Bi nt (r) due to lattice 
deformations, given by Eq. (19). 

The effect of strain in the low-energy electronic struc¬ 
ture of graphene can be accounted for by introduc¬ 
ing a scalar and a vector gauge potential in the Dirac 
equation 20,21,30 . The scalar term reads 

V w (r) = g[u xx (r)+u yy (r)], (18) 

while, for the K -valley and for an armchair crystallo¬ 
graphic orientation along the x-axis, A = (A x ,A y ) is 


FIG. 3. (color online) Extrinsic and intrinsic magnetic field 
correlations functions: (a) C e *t(r — r') = (B ex t(r)B ex t(r f )} 
for an external B\\ applied along the ^-direction and (b) 
Gint (j* — r') = (B[ n t(r)Bint(r')) due to strain in units of 
tfms /\ 6 (h/3 /ea) 2 , both corresponding to a ripple disordered 
surface h( r) characterized by a Gaussian correlation function. 

given by 

A x (r) = ^^[u xx (r) - u yy {r )], 

A y( r ) = - 2 ~ Ux v^ r ) > ( 19 ) 

where a = 1.42 A is the bond length between nearest 
neighboring carbon atoms, and e is the electron charge. 
Here g ~ 4 eV, k ~ 1/3, and /3 - —d\ogt/d\oga ~ 2 are 
dimensionless material dependent parameters 19,31 that 
characterize the coupling between the Dirac electrons and 
the lattice deformations, where t ~ 3 eV is the hopping 
integral between nearest-neighbor 7r-orbitals. 

In summary, for any given h(r) one can readily cal¬ 
culate the pseudo-magnetic field £?i nt = V x A. Since 
A z = 0 and neither A x nor A y depend on z, B- int = B- int z. 
Figure 2(c) shows the B[ nt corresponding to the random 
rippled surface h(r) of Fig. 2(a). Notice that the typical 
correlation length of B{ nt (r) is much shorter than that of 
h{r). 

Let us calculate (Hi nt (r)Hi nt (r / )) for a random Gaus- 
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sian correlated h(r), corresponding to Eq. (6) with 
F(x) = e~ x / 2A . To this end, we calculate the Fourier 
transform of the intrinsic pseudo magnetic field, namely 


with 


~ q'jMq'Mq - q'), ( 21 ) 


. h{3K 


ea 


We obtain 


h A 

1 L rrr 






where i and j label the Cartesian coordinates. 

We use Eqs. (20) and (21) to write the correlation 
function of R int in momentum space. The evaluation of 
(Bi n t(q)B[ nt (—q)} amounts to compute the correspond¬ 
ing (iq J -(< 7 )tq / ■,(— q)), that result in four-h correlation 
functions. This can be done exactly for Gaussian fluctu¬ 
ations and provides a good qualitative estimate for other 


-^int(^f) — ^ [Qy^xx (<?) T ‘^Qx'U'xy (<?) Qy^yy (O.)] •> (^0) Cases 


30 


16 + A 4 g 4 sin 2 3 0 
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( 22 ) 


where 6 is the angle between q and the x-direction. 

By Fourier transforming back to coordinate space, we arrive at 


Gnt(r 


nh = 


0 2 ^ 6 

20 & +9 f - 2 v sin23 “ 


2 / \ 2 
o~P /A 


(23) 


where a is the angle between p = r — r' and the x-axis. 33 

The correlation function C[ nt (r — r ') has 6 symmetry 
axes, reflecting the underlying graphene honeycomb lat¬ 
tice symmetry 19 . In other words, information about the 
graphene crystal structure survives disorder averaging. 
Figure 3(b) shows C{ n t(r — r f ) obtained from 10 5 numer¬ 
ical realizations of Gaussian correlated disorder for h(r). 
The numerical simulations serve as a helpful test to check 
our analytical results. As in the previous subsection, we 
verify an excellent agreement within the statistical pre¬ 
cision. 


III. DRUDE-BOLTZMANN CONDUCTIVITY 

In this Section we use the effective Dirac Hamiltonian 
of Eq. (1) to calculate the transport time and the Drude- 
Boltzmann conductivity of graphene monolayers in the 
presence of random gauge fields. 

High mobility graphene samples have typical electronic 
mean free paths of £ > 50 nm. Recalling 1 that the 
carrier d ensity is related to the Fermi wave number by 
kp = vMnI, one readily obtains that kpi >> 1 al¬ 
ready for a doping where \n\ ~ 10 11 cm -2 . This indi¬ 
cates that even for modest carrier densities a semiclas- 
sical transport description is justified. For \n\ > 10 11 
cm -2 the typical graphene conductivity in good sam¬ 
ples is much larger than e 2 //i, the order of magnitude 
of quantum contributions to the electronic transport, 
such as weak localization 2,26,27 and universal conduc¬ 


tance fluctuations 34 . In such situations, the Boltzmann 
approach is very successful in assessing the conductivity, 
as shown by direct comparison with numerical simula¬ 
tions using an atomistic basis 6,35 . As one approaches the 
charge neutrality point, and kpt < 1, the semiclassical 
method is no longer suited and one has to resort to more 
sophisticated approaches 2,8 . 

We now discuss how to add gauge field disorder in the 
Boltzmann approach. For long ranged disorder, some 
authors 36-38 argue that it can be advantageous to include 
the disorder potential in the classical Liouvillian evolu¬ 
tion, that is, to treat V(r) in the left-hand side of the 
Boltzmann equation. This approach is justified in the 
“classical” regime, where kp\ 1, that is, where the 
random fields with a characteristic length A vary slowly 
in the scale of kf}. In graphene samples with ripple 
sizes A of the order of few to ten nanometers 13-16 , the 
latter inequality holds for a carrier density \n\ 10 12 

cm -2 , that is much larger than the doping studied in 
most experiments 1 . 

In this study, we calculate transport times for both 
short and long ranged disorder by evaluating the corre¬ 
sponding Boltzmann collision integral (at the right-hand 
side of the equation). As mentioned in the introduction, 
for graphene on standard substrates (the case of interest 
here), the random gauge field contribution to the con¬ 
ductivity is not the dominant one. Hence, the electron 
mean free path due to ripples is larger than i and the ar¬ 
guments justifying the semiclassical approximation hold. 

The Boltzmann equation for graphene under a uniform 
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electric field E reads 1,39 


F de k , s df 0 
eE '~dP^ 


^ ^ ( 9 k,s 9 k',s')Wk',s<^k,s 

k' ,s' 


(24) 


where e k ^ s = svph\k\, fo is the Fermi distribution func¬ 
tion, Qk — fk ~ fo is deviation from equilibrium due to 
the electric field, and W k ^ s '^ k ,s is the transition rate 
from state (fe, s) to (fe', s), that we calculate using Fermi 
golden rule, namely 

O'Tr 

= -jr (|<fcV |V|fcs)| 2 ) 5(e k , s - e h >,s>), (25) 

where V is a generic long-ranged disorder potential 
parametrized by Eq. (2). The ^-function reflects the fact 
we are dealing with elastic processes and, hence, s = s '. 
In our model, the transition rates do not depend on s. 
Accordingly we drop this index whenever its omission 
does not introduce an ambiguity. 

The scattering processes we address are anisotropic. 
The calculation of the transport properties in this case is 
slightly different 40,41 than that of the standard isotropic 
case 39 . In this study, we adapt the nice method devel¬ 
oped by Tokura 42 - that is briefly described in what fol¬ 
lows - to calculate the transport times of massless Dirac 
electrons in graphene. 

In both situation of interest, the scattering potential 
correlation functions have at least one symmetry axis. 
For convenience, we choose the x-axis along a symmetry 
axis and define 42 

9k = (y eVk r ( 6 0 • ( 26 ) 


where r(0) is the relaxation time vector to be solved. We 
recall that 0 is the angle between k and the x-axis. Note 
that r depends explicitly on 0 and implicitly on \k\. 

The current density (spin and valley degeneracies in¬ 
cluded) is given by 


4 y- 

3 = ^ 2_^ eVkgk 

from which one obtains the conductivity tensor 

_ e 2 |g.F| f 27T , ft ( t x ( 0) cos9 t v {9) cos 0\ 

17 h 2 TT 2 J 0 \ t x (0) sin# Ty(0) sin# ) ’ 


(27) 


(28) 


where Sp is the Fermi energy, measured with respect to 
the charge neutrality point energy. For the sake of sim¬ 
plicity, in Eq. (28) we have taken the zero-temperature 
limit, namely, —dfo/de = S(e — £f). 

By substituting the ansatz (26) in the Boltzmann equa¬ 
tion, Eq. (24), one obtains an integral equation for t( 0), 
namely 

/»27T 

cos6 = / d0'[T x (0) - t x (9')]V9{0,9') (29) 

Jo 

sin 0= / d0'[Ty(9)-T y (e')}W(0,e'), (30) 

Jo 


where 9' is the angle between k' and the x-axis. and 


W(M') 


— f 

(27T) 2 Jo 

A\sf\ 

2pv 2 F h? 


dk k Wk',s<—k,s 

(| (k’,8\V\k,8 )\ 2 >, 


(31) 


where, due to the zero-temperature limit, k = k! = kp. 

The matrix element (fe', s\V\k, s) depends on q = k — 
k' and p = 7r/2 + (0 + O')/ 2, the angle between q and 
the x-axis. Hence, W(0, 0') is better cast as W(q^p). 
We use the standard notation £ = \0 — 0'\ and = 
2 kp sin(£/2). 

By expressing r x and r y in terms of a Fourier series, 
one transforms Eqs. (29) and (30) into an (infinite) set 
of algebraic equations. Using the x-axis symmetry and 
that W(0,0') = W(0',0) 43 , one shows that 42 

oo 

Tx(0) = E 4 n) cos[(2n - 1)9] (32) 

n= 1 
oo 

Ty(9) = J2^ n) M(2n-m. (33) 

n= 1 

By inserting the above relations in Eq. (28), one con¬ 
cludes that the conductivity tensor is diagonal, with 


^ XX 


e2 M (i) 

H 2 TT x 



(34) 


that supports the interpretation of r^ as a transport 
time vector. 

The symmetry W(0, 0') = W(0',0) implies that 

W(qp p) = yV(qp p + 7r) = W(qp p — 7r). In turn 


w(g c , <p) = E w «(9c) cos(2n<p), (35) 

n=0 


with an obvious inversion relation. 

By using the Fourier expansions for r(0) and W (g<y p), 
Eqs. (29) and (30) can be cast in matrix form 42 

oo oo 

*M = E M »7n4 n) and 5 lA = J2 M JnP\ (36) 

n= 1 n=1 

where the matrix elements of M ± are 42 

(_ 1 y—n 

n — [(1 T ^Z,n)^|/-n|,n+«-l ^ ^n+l — 1 s |Z—n|] i 

(37) 

with 

Jn,m= / dCW n (g c )[cos«) - cos(mC)]. (38) 
Jo 

Finally, by inverting in Eq. (36), one writes the vec¬ 
tor transport time components as 

4 1} = [(M“) _1 ]n and rfjp = [(M + ) _1 ]n. (39) 
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Note that for isotropic scattering, all W n with n > 
0 are zero and M ± is diagonal, with elements Kij = 
Jo 5 2 Z-i- Hence, the vector transport time components 
coincide, r^ ^ = t^ 1 ), and read 

i r 2n 

= Jo,i = j ^ ^CW 0 (gc )( 1 - cosC), (40) 

which is the standard expression for the transport time 
in isotropic systems. 


A. Effect of an in-plane magnetic field 


Let us now calculate the effect of an external parallel 
magnetic field on the conductivity. From Eq. (13) we 
write the effective disorder potential for the A-valley as 


V ext (r) = v F ecjyAy{r) = -v F eB\\h(r)a y . (41) 

We recall that h(r) varies slowly in the scale of the lattice 
spacing and, hence, V ext (r) is long-ranged and does not 
mix valleys. 

The momentum relaxation rate Wk'^k reads 


MWfc = 5(k-k , ) 2 ^y^ Bf l (42) 


where 


C h (q) = Jdr e^ r (h(0)h(r)) (43) 


is the form factor of the height-height correlation func¬ 
tion. Here (• • •) indicates disorder average. 

From Eq. (31) we obtain 

>v(<?,y>) = r^- C h {q)(l + cos2y>), (44) 

that has only 2 non-zero Fourier components, namely, 


W„(<?) 


(eg||) 2 M 

47 rh 3 


C h (q ), 


for n = 0,1 


(45) 


while W n (q) = 0 for n > 2 . 

In this case, the M ± matrix is tridiagonal and reads 42 


M ± 


((lT^o,i 0 • • 

— Jo, 3 

0 «/o,5 

V 


(46) 


The inverse transport time components are given by 

13 11 

~(iy = 2^°4 _ ^3 and -^y = - Jo,i — F 3 , (47) 

Tc Ty 

where T 3 is isotropic and determined by the continued 
fraction relation 


r 


m 


(Jl,m-l ) 2 

4(J 0 ,m r m+2 ) 


(48) 


In practice, we compute T 3 by assuming that T^ = 0 and 
subsequent iteration of (48). The choice of m determines 
the precision of the calculation: The larger m, the more 
accurate is T 3 . 

Assuming that T 3 Jo,i leads to an interesting result, 
that is 

“ ii) = (eB sl» F ' l ‘Kd-cosOCM- (49) 


In this limit Ty^ /r^ = 3. In other words, for T 3 <C Jo,i 
the corrections to the conductivity due to B\\ lead to 
A G yy — 3A a xx , regardless of the dependence of correla¬ 
tion function Ch(q) on q. 

For r 3 = 0 and Ch(q ) = 27tA 2 h 2 Yns e~ x2q2 / 2 , we write 
Tx]y in closed analytical form, namely 


r (!) 


3r £ 


(i) 


(eB ||) 2 | 5 F | 
4 h 3 


(A h rms ) 


xe- A ^[/ 0 (A 2 4)-/ 1 (A 2 4)], (50) 


where Iq and I\ are mo dified Bessel functions of the first 
kind. We use 1 k F = \fr\n\ to express the conductivity 
in terms of the charge carrier density n. We conclude 
that the correction to the conductivity due to an in-plane 
magnetic field depends quadratically on h rms B y and has 
a non-trivial dependence on A 2 |n|. 

In the high doping limit of A Ini 1 / 2 1, Eq. (50) gives 
the resistivity contribution of an in-plane magnetic field 
in a rippled graphene sheet 


J^Pyy — 3A p xx 


1 drills 

2y/2 h A 



(51) 


in agreement with Ref. 9. 

Figure 4 shows the resistivity A p yy versus the car¬ 
rier density n (due to particle hole symmetry, we only 
show n > 0) for m = 3 and m —)> 00 . The optimal m 
value to obtain convergence depends on A 2 n. The inset 
shows A pyy for A 2 n values outside the validity range of 
the asymptotic expansion. As discussed in the next sec¬ 
tion, the A 2 n range displayed in the inset corresponds 
to the typical experimental situation. We find that the 
|n | -3 / 2 scaling predicted by the asymptotic expansion 
(51) is only a rough approximation. 

In general, T 3 is a non-vanishing correction to the 
transport time components, hence /t x ^ 7 ^ 3. How¬ 
ever, for Gaussian ripple height correlations, the ratio 
Ty^/r^ is a function only of A k F . 

Figure 5 shows Ty 1 ^ /t x ^ versus \ 2 n. It illustrates the 
importance of T 3 in the calculation of the conductivity 
corrections. We find that by increasing the carrier con¬ 
centration the anisotropic conductivity is considerably 
favored. 

The transport properties predicted by Eq. (50) are only 
slightly modified for the case of exponential ripple heigh 
correlations, (h(0)h(r)) = h 2 ms e -r / A : For AItt ) 1 / 2 >> 1, 
the resistivity tensor given by Eq. (51) is multiplied 44 by 
a prefactor of order of unity times log(A 2 |n|). 
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due to strain that has been reported in the literature 30 
is not sufficient for the analysis we propose. 

As in Sec. IIB, we calculate ( Uij(r)ui/jt(r ')) by as¬ 
suming Gaussian correlated ripple height fluctuations, 
Ch(q) = 27rA 2 h^ ms e~ x2q2 ^ 2 . After some lengthy but 
straightforward algebra, we obtain 


<|<fc' s |K lt |fcs}| 2 ) 


Vph 2 /3 2 K 2 7rhf mi 

32 a 2 AA4 

3 


(53) 


x 116 + A A q A cos 
Using Eq. (31) we arrive at 


2 ^ + 0 ') 


a -AV/4 


FIG. 4. (color online) Resistivity A p yy in units of po = 
(7rA/i rm sB||) 2 /2^ due to B\\ as a function of A 2 n. Inset: The 
same as in the main plot in log-log scale to illustrate the de¬ 
pendence of A p xx on \n\. 


W(q,<p)=W int 


A 4 4 

16 H- ^-(1 + cos 6 <^) 


„-AV/4 


, (54) 


where 



FIG. 5. (color online) Ratio r y ^ /r^ as a function of tt\ 2 n. 
For 7rA 2 n < 5.0, the limit m oo is attained within 10 -4 
accuracy for m — 13. 


B. Effect of strain fields 

Let us now consider the effect of the pseudo-magnetic 
fields due to strain, J?i nt , in the conductivity of monolayer 
graphene sheets. In contrast to the mechanism discussed 
above, B{ nt (r) is solely determined by h(r) and the mate¬ 
rial properties. Hence, it intrinsic to any graphene sample 
with disordered ripples. 

We use Eq. (19) to calculate the vector transport time 
r for the intrinsic effective vector potential due to strain. 
For the K- valley Unt reads 

Vint(r) = hv F — | [u xx (r) - u yy (r )] a x - 2 w xy (r)< 7 3/ |. 

(52) 

In contrast with the previous subsection, here it is dif¬ 
ficult to make quantitative progress without assuming a 
specific form for the ripple height correlation function. 
The qualitative behavior of the conductivity corrections 


= /? 2 /t 2 |e F | fr 4 ms 

64:a 2 h A 2 


(55) 


(The notation is the same as that of the previous section.) 
The only non-zero Fourier components of W(g, cp) are 

Wo (q) = Wint (f 6 + ) e _A2<?2/4 and 

W 3 ( 9 )=W int lh!e- AV / 4 . (56) 

In this case, the M ± matrix reads 42 


M ± = 


Jo.l 0 t|J 3 ,2 -\Jza “X 

0 Jo,3 T \ J 3 ,o 0 0 

T5J3.2 0 Jo,5 0 

— 2 ^ 3,4 0 0 Jq,7 


V 


7 


(57) 

and the inverse transport time components are given by 


■JIJ - J 0d _ 

x/y 


(J 3 , 2 ) 2 (J 3 ,4) 2 

4 Jo 5 4 Jo 7 




(58) 


Since t^ = r y ^, the conductivity corrections due to the 
strain field are isotropic. This result seems to be at odds 
with the fact that the pseudo-magnetic field autocorre¬ 
lation function (Bi nt (r)Bi nt (r')) clearly shows an hexag¬ 
onal symmetry, as illustrated by Fig. 3b. As shown by 
Tokura 42 , using general arguments, this is a false para¬ 
dox: The conductivity tensor becomes anisotropic only 
for scattering processes characterized by a single symme¬ 
try axis, like in the B\\ case, analyzed in the previous 
subsection. 

Assuming that Joq vastly dominates the sum in 
Eq. (58), we obtain 
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r (1) 

x/y 


= Wi„ t 7re 


-X 2 k 2 F /2 


(32 


■ 8X 2 kp 


16X 4 kp)Io(X 2 kp/2) - (64 


24 A 2 k 2 p 


16X 4 kp)Ii(X 2 kp/2) 


(59) 


where / (J and I\ are modified Bessel functions of the first 
kind. 

In the limit of A|n | _1 1, we obtain the correction to 

the resistivity 


An -An ss A 2 ^ 2 ^ 2 ^rms ^-3 11 - 3/2 (601 

Apx X -Apyy ~ e 2 32?r a2(j2 A \n\ . (bUj 

The above asymptotic leading order expansion for A p 
helps us to develop some insight on the relevant param¬ 
eters. However, since the situation of kp\ 1 is hardly 
met in the current experimental situations of interest, it 
is necessary to numerically calculate the inverse transport 
times, as we do in the next subsection. As expected the 
strain corrections to the conductivity depend on material 
parameters, and are a non-trivial function of A, h rms , and 
\n\. Since these corrections are small compared to other 
disorder effects, they are difficult to be noticed in stan¬ 
dard transport experiments. This situation changes if 
we consider the combined effect of intrinsic and extrinsic 
random magnetic fields, as we discuss in the next. 


C. Combined effect of extrinsic and intrinsic 
random magnetic fields 


tern transport properties are dominated by other scatter¬ 
ing processes, with a corresponding (isotropic) transport 
time r s . 

It is customary to use Matthiessen’s rule when dealing 
with systems characterized by different competing relax¬ 
ation time mechanisms. In our case, Matthiessen’s rule 
translates into adding the inverse transport times given 
by Eqs. (47) and (58), namely 



2 zb 1 


(4D 2 

441 


■ Text -pext I Tint 
" J 0,l — 1 3 “I" ^0,1 — 


(4D 2 

44 n 7 

(61) 


This naive approach was shown to be inaccurate when 
dealing with anisotropic potentials 42 . 


We analyze the combined effect of intrinsic and extrin¬ 
sic random magnetic fields by considering an effective 
M-matrix given by 


M± t = M± t + M± t , (62) 

where M^ t and M± t are given by Eqs. (46) and (57) 
respectively. 

The r components are 


We conclude this section by analyzing the combined ef¬ 
fect of both previously discussed sources of random mag¬ 
netic field disorder. As before, we assume that the sys- 


4 1} = K M tot ) 1 ]ii and 4 1} = [( M tot) 1 ]ii, (63) 

that, for m = 5, explicitly read 


x/y 


Jo, 


int 


(J$$ + 41T 41/2) (415 + 41) - (41)74 


( Text , Tint \ ( Text \ 2 
A WO, 5 ~r t '0,5/ W 1,2 ) 


4 ‘ 


Text Text 

3 2 2 J 1A 


- 4 ( 4 ( 2 ) (4 


ext 

3 


Tint 

J 0,3 


T 241) 


(64) 


This result is clearly different from Eq. (61), since it 
mixes intrinsic and extrinsic effects. Let us now dis¬ 
cuss the dependence of r on B\\, A, h rms , and n. For that 
purpose we numerically invert the matrix M t Q t , at order 
fn ~ 30 — 50 to guarantees an accuracy of 10 -5 for the 
analyzed parameter range. 

The resistivity corrections obtained from the 
Matthiessen’s rule, Eq. (61) depend quadratically 
on H||, in line with the experiment 9 . The full M-matrix 
analysis does not guarantee this simple dependence. 
In Fig. 6 we plot the resistivity correction A p yy cal¬ 
culated using the full M-matrix and compare it with 
the one obtained from the Matthiessen rule, given by 
Eq. (61), for realistic values of h rms , A, and n. The full 


M-matrix calculation (indicated as “exact”) shows an 
overall higher resistivity than that obtained from the 
Matthiessen rule. It depends linearly on E > 2 for large 
B\\ and deviates from this dependence only when B\\ 
becomes small. 

In contrast, the dependence of Ap with h rms , A, and 
\n\ is not trivial. For h rms and A values taken close to 
the ones reported by topography experiments 13-16 , a nu¬ 
merical study using the full M-matrix approach gives 
A pyy oc A -a with a ~ 3 • • • 4, A p yy oc hf ms with /3 ~ 3, 
and A p yy oc \ri\~ 1 with 7 « 2 . In summary, Ap is very 
sensitive on small variations of h rms and A. 

In Fig. 7 we compare Eqs. (47), (61), and (64) to gain 
insight on how the strain mechanism affects the ratio 
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FIG. 6. (color online) Resistivity correction A p yy in units of 
po = (7r\h rms B\\) 2 /2h as a function of the in-plane magnetic 
field B\\ for h r ms — 0.2 nm, A = 30 nm, and n — 10 12 cm -2 . 
Inset: The same as in the main figure for h r ms = 0.4 nm and 
A = 10 nm. 

Ty 1 ^ /t£~\ We find that the strain fields contribute to 
a strong suppression of the anisotropy in the transport 
time due to a strong B \\. However, for realistic parameter 
values the anisotropy is still very large and of the order 
of Ty 1 ^/t^ ~ 10 for \n\ ~ 10 12 cm -2 . 



FIG. 7. (color online) Anisotropy r^/ri^ as a function 
of the carrier concentration n using different approximation 
schemes, for /i rms — 0.2 nm, A = 30 nm, and B\\ — 8 T. The 
red and green lines stand for the contribution of B\\ with¬ 
out accounting for strain fields. The blue line represents the 
contributions of both external and strain fields using the con¬ 
ventional Matthiessen’s rule. The black line stands for the 
combined effect of intrinsic and extrinsic fields obtained for 
the full M -matrix analysis. 

In order to further compare our results with the 
experiment 9 , let us introduce the magnetoresistance 
A p = E ■ j/f, where j x = j cos£, j y = j sin£, and 
£ is the angle between B\\ and j. Using the relation 
Ei = pijJj one writes 45 

A p{£) = A p xx cos 2 £ + A pyy sin 2 £. (65) 


Ref. 9 reports Ap(70°)/Ap(20°) « 0.13 - 0.26. Using A 
and h rms values obtained from AFM measurements, we 
obtain T^/ri^ ~ 10 for n — 10 12 cm -2 . This ratio leads 
to Ap(70°)/Ap(20°) ~ 0.2 in good agreement with the 
experiment 9 . 


IV. CONCLUSIONS AND OUTLOOK 

In this paper we studied the effect of random magnetic 
fields on the transport properties of a rippled graphene 
flake. We used the Boltzmann equation, adapted to the 
case of anisotropic disorder 42 , to address the case of an 
external magnetic field applied in-plane, the effect of in¬ 
trinsic strain fields caused by the graphene corrugation, 
as well as the combination of both. 

We find that an external in-plane magnetic field B\\ 
gives rise to very anisotropic conductivity corrections. 
By neglecting the effect of strain fields and using a 
parametrization of the ripple disorder that is consistent 
with experiments, we find conductivity corrections that 
scale with B jj and |n| -2 , consistent with Ref. 9. In con¬ 
trast, we obtain r^/ri^ ratios as large as 20 • • • 30. 

In the absence of an external magnetic field, random 
gauge fields due to ripple disorder give a isotropic contri¬ 
bution to the electron momentum relaxation in graphene, 
in line with the order of magnitude estimate presented in 
Ref. 30. We find, however, that ripple disorder cannot 
be neglected in the analysis of the conductivity in the 
presence of a large B \\. We also conclude that, due to the 
anisotropic nature of the problem, Matthiessen’s rule is 
not accurate to address both intrinsic and extrinsic ran¬ 
dom fields at the same footing. For that purpose we have 
to invert the total M-matrix. 

This approach allows us to successfully describe the 
corrections to the Drude conductivity reported in the 
experiment 9 using typical A and h rms parameters taken 
from the AMF literature. In addition, we also obtain a 
suppression of the resistivity anisotropy (with respect to 
the case where strain is neglected) that consistent with 
Ref. 9. We believe that the anisotropic nature of the ran¬ 
dom magnetic field disorder also significantly changes the 
quantum corrections to the conductivity with respect to 
the isotropic result 28 and deserves further investigation. 

In summary, our results suggest that the investigation 
of anisotropy corrections to the Drude conductivity can 
be a new and insightful path to experimentally quantify 
effects of random pseudo-magnetic fields due to strain. 
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